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Abstract 

We use density-functional theory, of the fundamental-measure type, to study the relative stability 
of the biaxial nematic phase, with respect to non-uniform phases such as smectic and columnar, in 
fluids made of hard board-like particles with sizes cJi > cj2 > (J3. A restricted-orientation (Zwanzig) 
approximation is adopted. Varying the ratio ki = ai/(J2 while keeping K2 = o'2/cj3, we predict 
phase diagrams for various values of K2 which include all the uniform phases: isotropic, uniaxial 
rod- and plate-like nematics, and biaxial nematic. In addition, spinodal instabilities of the uniform 
phases with respect to fluctuations of the smectic, columnar and plastic-solid type, are obtained. 
In agreement with recent experiments, we find that the biaxial nematic phase begins to be stable 
for K2 — 2.5. Also, as predicted by previous theories and simulations on biaxial hard particles, 
we obtain a region of biaxility centred on ki ~ K2 which widens as K2 increases. For K2 ^ 5 the 
region K2 ~ ki of the packing-fraction vs. ki phase diagrams exhibits interesting topologies which 
change qualitatively with K2. We have found that an increasing biaxial shape anisotropy favours 
the formation of the biaxial nematic phase. Our study is the first to apply FMT theory to biaxial 
particles and, therefore, it goes beyond the second-order virial approximation. Our prediction that 
the phase diagram must be asymmetric is a genuine result of the present approach, which is not 
accounted for by previous studies based on second-order theories. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



The existence of biaxial nematic (Nb) phases is a problem of great fundamental and 
practical interest, since these phases could be used in fast electro-optical devices hj. After 



lOj and experimental observation. The Nb phase was detected in a 



yotropic fluid [ll| and, 

more recently, in liquid crystals made of bent-core organic molecules |12|-|l5l|: however, in the 
latter case controversial issues about the correct identification of the Nb phase still remain l|, 
1^ 17 1 . Biaxial phases have been predicted 18 1 and observed in computer simulations [lol of 



mixtures of rod-like and plate 



and demixing is expected 



20- 



ike particles. However, strong competition between biaxiality 



24| . and demixing could preempt biaxiality altogether. 



Recently, a promising approach based on colloidal dispersions of mineral particles has been 
proposed. Goethite particles with a board-like shape of sizes ai > a2 > as and exhibiting 
short-ranged repulsive interactions were synthesised and observed to form colloidal liquid- 
crystalline phases of nematic, smectic and columnar type 
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nematic phase has been observed 



28 



27|. Interestingly, a biaxial 



when (Ji/a2 ~ o'2/o'3 = 3. In this case the Nb phase 
seems to be more stable than the competing smectic and columnar phases. Indeed, theory 
and simulation predicted that the biaxial nematic phase will be formed when ai/a2 ~ 0-2/0^3 

The existence of biaxial phases in these colloidal particles opens a new avenue for theo- 
retical research. Their short-ranged repulsive interactions and well-defined shape allows for 
the application of sophisticated density- functional theories for hard particles. In particular, 
these studies may help to understand why the biaxial phase is so hard to stabilise and if 
competition with other more highly ordered phases, such as smectic and^columnar, could be 
at work. In this connection, the early work by Somoza and Tarazona 
weighted-density functional theory applied to hard oblique cylinders, found stable smectic 
phases and only a very limited region of biaxial- nematic stability. In more recent work, 
Taylor and Herzfeld SjJ] considered a fluid of hard spheroplatelets (which are similar to 
board-shaped particles but with rounded corners) and examined its phase behaviour using 
the scaled-particle theory of Cotter 32|, which incorporates the excluded volume associated 
with two particles, combined with a free-volume theory. The predicted phase diagram again 
showed that biaxial-nematic stability was almost entirely suppressed by the smectic phase. 
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30|, who used a 
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A recent simulation work 



33j analysed the full phase diagram of board-like particles with 



square cross section, ai = and 0.125 < o-ijax < 0.5. The first condition is not compati- 
ble with the formation of the the biaxial-nematic phase. To our knowledge, no simulation 
studies have been undertaken so far concerning more suitable size parameters. 

Here we report on a density-functional theory where the phase stability of a fluid of 
hard board-like particles is explored. Our aim is to examine the phase behaviour of the 
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focusing on the 



goethite-based particles used in the experiments of van del Pol et al. 
formation of the Nb phase and their interplay with the uniaxial nematic (Nu) phases and 
the smectic (S) and columnar (C) phases. The approach is based on a fundamental-measure 
theory (FMT) for hard parallelepipeds 3J] in the so-called Zwanzig or restricted-orientation 
approximation. Using this theory, we predict islands of stability for the Nb as a function 
of the size parameters. These islands are not based on rigorous boundaries, as the effect 
of the S and C phases is estimated only via bifurcation analyses (not via costly free-energy 
minimisations). However, our work should provide trends as to which particles sizes are 
more optimal to stabilise the Nb phase. As a general rule, the particle aspect ratios should 
be larger than those explored in the experiments. Our work extends and improves a previous 



theoretical investigation 



3l| which considered very limited particle aspect ratios. Also our 



study, which is the first to adopt the FMT approach for biaxial particles and, in particular, 
for hard board-like particles, may have some relevance in display applications based on 
fast-switching particles that reorient by rotating their short axes instead of the long ones, 
something that can be achieved with biaxial nematic phases. Finally, polydispersity in the 
sizes may play a key role when comparing our results with the experiments, but was not 
taken into account in the present work. 

In Section [TTl we define the particle model, and summarise the FMT theory and order 
parameters used to characterise uniaxial and biaxial order. Also, we write the bifurcation 
equations solved to obtain the spinodal lines that define instability against non-uniform types 
of ordering. Section [IHl presents the results, and in Section [IV] we give some conclusions. 
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II. MODEL 



The system consists of a collection of biaxial hard parallelepipeds with edge lengths (jj 
(z = 1, 2, 3) satisfying ai> a2> 0-3 (Fig. [T]). We define the size ratios 

«i = — , 1^2 = — , (1) 

as the ratio of length-to-width and width-to-thickness, respectively. Now the fluid of iden- 
tical particles is mapped onto an equivalent fluid mixture of six species as follows. Within 
the Zwanzig approximation, where particle orientations are restricted to the three orthogo- 
nal Cartesian axes, particles oriented in the same way are considered to belong to the same 
species, and six unequivalent species result. Then, the density profiles p^yijc) will be the fun- 
damental quantities that characterise the equilibrium properties of the fluid; the subindices 
/i and V represent the directions parallel to the first (length) and second (width) main axes 
of the particles, respectively. 
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FIG. 1: Schematic of board- like particle with definition of the three sizes: ai (length), (72 (width) 
and (T3 (thickness). A depiction of the six possible particle orientations is included. Labels iiv 
(with = {x, y, z} and ^ v) indicate orientation with longest axis parallel to ^ and second 

axis along u. 



According to FMT 3J| , the excess part of the free-energy density in thermal-energy units 
kT = is given by 

(2) 



^ 1/1 X ni ■ n2 , n2xn2yn2z 



l-ris (1-^3)2' 
where the weighted densities na{r) are defined through convolutions: 



(3) 



with a = 0,1,2 or 3. The weights a;^y''(r), defined in 3^ for uniaxial parallelepipeds, are 



extended in the present work to the general case of biaxial particles as 



^ 3 3 

^U^^uix.). u;(lHr) = nV(^0, (4) 



i=l i=l 



5f,u{,x)5^,y{y)e^y{z)) , (5) 



2 

0,.{x)9^M^^u{z))^ (6) 



where we have defined 



5i,u{xi) = 6 (^-^ - \xi\j , (7) 

O^^uix,) = e(^^-|x,|^, (8) 

with 6{x) and G(x) the Dirac-delta and Heviside functions respectively, while o"^^ = o"3 + 
(cTi — o'3)6i^ + ((72 — o'3)6ii, denotes the length of the biaxial particle with orientations fi and 
u along the z-th direction. We have used the notation xi = x, X2 = y and xs = z. 

The total free-energy density is obtained by adding the ideal free-energy contribution, 

I3F = ! rfr[<l>id(r) + $(r)], (9) 
Jv 

$M(r) = J]p^.(r)(lnp^,(r)-1), (10) 

where V is the fiuid volume. The equilibrium state of the fiuid follows by minimisation of 
PJ^ /V with respect to the six independent densities p^u{^) at fixed total density per unit 
cell (the latter being one-, two- or three-dimensional for the smectic, columnar or crystalline 
phases, respectively). 

A. Uniform phases: isotropic, uniaxial and biaxial nematics 

In the uniform phase the free-energy densities can be written as 

$id = p I Inp - 1 ^X/,i.lnx^i. I (11) 

$ = -eoln(l-e3) + f4 + |^^, (12) 

1-6 (l-Cs)^ 



with = P^^y^nuW^/w ■ The total density and packing fraction are = P and ,^3 = 77, 
respectively, while x^i, are the molar fractions. The coefficients wjS) are defined as w^f^ = 
J drools) (r), with w^^J = 1, wj^J = 0-10-20-3, and 



w 

fj,i/ 



il^ = fc, (13) 



.(2) 



w;.;' = • (14) 

The fiuid pressure can be calculated from 
while the chemical potential is 

/3/x = p-i(/3p+$id + $). (16) 

To find the equilibrium uniform phases we fix the packing fraction rj and minimize the total 
free-energy density $id + $ with respect to the molar fractions {x^i/}, with the obvious 
constraint "^^^x^^, = 1. 

Let us now discuss how the orientational order can be defined quantitatively. In the 
present model, the usual order parameter tensor 

Qap = ^ (3(m>P - Sa^) , (17) 

(with the projection of the unit vector of the long axis of particle j on the a Cartesian 
axis, the average taken over all particles) becomes diagonal, 

Qal3 ~ 2 ( XfJ.uSlJ.aSlJ.l3 ~ j 

= I (^Sj^x^^-lj 6^1,. (18) 

Assuming that the first nematic director is parallel to the z axis while the second lies along 
the X axis, our uniaxial and biaxial order parameters are, respectively 



Q = Qzz = ^ Xzu - 1^ 



S — Qxx Qyy r) [ ^ ] -^xu ^ ] Xyi, | . (19) 



The order parameter S measures the difference in the fraction of particles pointing along the 
X and y axes. If there are no particles in the xy plane, then 6 = 0. Therefore, the order 
parameter matrix has the usual diagonal form 



/ Q-6 



Q 



\ 










Q + S 









QJ 



(20) 



However, the biaxial order parameter S, which is zero for perfect unaxial order and non-zero 
when there is some amount of biaxiality, has the problem that, for perfect biaxial alignment, 
S = 0. A better tensor quantity to characterize biaxiality involves the other two particles 
axes: 



(21) 



where and are the projections of the second and third axes of the j-th particle on the 
a Cartesian axis, respectively. For our system we again find a diagonal tensor. 




X 



6. 



al3- 



(22) 



These magnitudes can be writen as a function of new order parameters S and A as 



(5 -A) 



A 



(^ + A) 



yy 



A. 



S. 



(23) 



In terms of the mole fractions these order parameters can be explicitly written as 



S 
A 



•'^xy ~l~ -^yz •^yx ) ) 



1 . 

2 \-^xz 

1 

{•^yx -^xy -^xz -^yz ~^ '^{p^zx 



zy) 



(24) 
(25) 



Note that the first and second nematic directors have been chosen to be parallel to the z 
and X axes, respectively. While S measures the extent of biaxial ordering along the x and y 
axes, A is the global biaxial parameter which is zero for the uniaxial nematic phase (since 



•^xyy -^xz 



Xyz and 



uniaxial nematic symmetry implies Xyx 
perfect biaxial ordering (for which x^i, = 0, V/x, u z,x, while x 



X 



zy)i 



and unity for 



!)• 
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B. Non-uniform phases 



The existence of spatially non-uniform phases, such as the smectic or columnar phases, 
has been assessed via bifurcation analysis from a uniform (either isotropic or nematic) phase. 
The relevant response function of the system is the direct correlation function, c^j,_T-6i(r — r'). 
We define the 6x6 matrix S with elements 



where a = /i + z/ — 2 if /i<z/ and a = /i + z/ + lif/i>z/ while /3 = r + ^ — 2ifr<^ and 
P = r + 9 + lifr>6. The Fourier transforms of the direct correlation functions. 



with q the wave vector, can be calculated from the Fourier transforms (q) of the corre- 
spondings weights a;^"^(r). To calculate the spinodal or bifurcation curves, we need to solve 
the equation 



where Dg = det (S). From here we find, by iteration, the values of packing fraction rj and 
wave vector q at the spinodal instability. Note that, during the iterative process, for each 
value of rj one needs to find the molar fractions x^^ via minimization of the free-energy 
density. 

III. RESULTS 

We have applied our FMT approach to study a fluid of board-shaped particles with fixed 
values of the width-to-thickness ratio, i.e. K2- Phase diagrams in the plane tj-Ki^ were 
produced in each case. 

Our first result is Fig. [2|^a), which depicts the phase diagram for the case K2 = 2.5. We 
identify one isotropic phase (I), with Q = A = 0, and two uniaxial nematic phases, N^^ 
and Nu"*, with Q ^ and A = 0. These two phases can be distinguished by the direction 

(t) 

the uniaxial director takes with respect to the sides of the particles. In , the director 
points along the long particle side ('rod-like' nematic), while in N^^ the director lies along the 



(26) 




(27) 



D,(r/,q) =0, VZ},(r/,q) = 0, 



(28) 
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FIG. 2: (a) Phase diagram in the plane rj — k~[^ for hard parallelepipeds with K2 = 2.5. (b) A zoom 
of the region about ~ 0.4. Labels indicate regions of stability of the different phases, or type of 
ordering with respect to which bifmxation cm'ves are computed (see text for key to labels). Solid 
and dashed curves indicate first and second order phase transition, respectively. Dotted curves are 
bifurcation lines to the corresponding non-uniform phase from the stable bulk uniform phases. In 
(b) the shaded region corresponds to the estimated maximum region of stability of the Nb phase. 



short side ('plate- like' nematic). In between these two phases a region of biaxial nematic, Nb, 
appears as a wedge centred at ki ~ K2, as predicted previously 3, ^, 6|. This region opens 
up as density is increased. Concerning the I-Nu transitions, they are always of first order, 
except at some particular point in the 77 — ^1 plane, where the density gap is exactly zero. 
This is the point where the continuous N^^'-Nb and Nb-N^'' transition lines meet (Landau 
point). These transitions are always continuous, regardless of the value of ki. The topology 
of the phase diagram, as far as isotropic and nematic phases are concerned, is similar to 



that obtained previously by Camp and Allen 6|] 
slightly different from that of Taylor and Herzfeld 



or hard biaxial ellipsoids. However, it is 



3l| for hard spheroplatelets since, in that 



case, the Nb region does not extend down to the I-Nu transition boundary. The restricted- 



orientation approximation was invoked in 



3l| to explain the discrepancy with simulation. 



However, the fact that we are also using this approximation points to a different explanation. 
The different particle geometries (a spheroplatelet is generated by a sphere whose centre is 
constrained to a rectangle, and therefore has rounded instead of sharp corners) may be one 
possibility. 
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One interesting feature of the -Nb transition line is that it tends to a constant packing 
fraction as — !■ oo. As one approaches this hmit, the value of the uniaxial order parameter 
Q tends to unity, and the transition exclusively involves a two-dimensional ordering of the 
intermediate (secondary axis) particle side. Also, in the neighbourhood of the value ki = K2, 
the biaxial phase becomes reentrant, as can be inferred from the negative slope of the Nb- 
visible in the zoom presented in Fig. which later becomes positive. 

Overall the phase diagram is very asymmetric with respect to the condition ki = K2, in 
contrast with models based on the second virial coefficient. This coefficient has rod-plate 
symmetry, and it is the inclusion (implicit in the present model) of third and higher virial 
coefficients that breaks the symmetry, producing an asymmetrical phase diagram. 



0.5| . 1 . TT". — I — ' — I — ' — I 0-182 




0.2035 



^1 ^1 

FIG. 3: (a) Phase diagram in the plane rj — ki for hard biaxial parallelepipeds with K2 = 5. (b) A 
zoom of the region about ~ 0.2. Labels as in Fig. [2j 

Let us turn to the spinodal lines. These are associated with instabilities of the stable 
bulk phase with respect to spatially non-uniform types of ordering. In Fig. [21 different 
bifurcation curves have been represented, corresponding to PS (plastic solid, with three- 
dimensional spatial order but no orientational order), Sy'' and Sy'* (uniaxial rod- and plate- 
like smectic order, respectively), Sb (biaxial smectic order) and C^^ (uniaxial columnar order 
with the long particle axes parallel to the column axes). In the region ki > 1 the I fluid is 
unstable with respect to the PS phase. This spinodal curve crosses the I-Ny'' coexistence 
region and connects to the S^ ^ spinodal curve, which flnally joins the stability region of the 
Nb phase. The region where the ^ phase is stable is reduced to a very small triangular 
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region. Inside the Nb region, we found a short bifurcation hne associated with Sb ordering 
[see enlarged area in Fig. [2](b)]. This hne touches the left boundary of the Nb region, and 
smoothly continues into the uniaxial region as a bifurcation line associated with columnar- 
type fluctuations, C^'' in Fig. [2] 35|, which in turn is connected to the S^^ bifurcation curve. 
These calculations indicate that, for the case K2 = 2.5, the region in which the Nb phase 



is stable is exceedingly sma. 
by Taylor and Herzfeld 



1. This coincides with the conclusion drawn for spheroplatelets 
3l|, who used a version of free- volume theory to study absolute 
stability of smectic, columnar and crystal phases. They concluded that the biaxial nematic 
phase is almost completely preempted by the uniaxial and biaxial smectic phases. 

The case k,2 = 2.5 was chosen to be close to that considered by Taylor and Herzfeld 



3l|. If a is the radius of the sphere and c x b are the sizes of the rectangle that define 



the spheroplatelet, with c > b > a, then the true sizes of the particle used in Ref. [3l| are 
(c + 2a) X (6 + 2a) x 2a, and we can define length-to-width and width-to-thickness ratios 
for this particle as Ki = (1 + c/2a) /(I + b/2a) and K2 = 1 + b/2a. In their work, Taylor and 
Herzfeld fixed the ratio c/2a = 5 and varied b/2a. This is equivalent to fixing the product 
K1K2 = 6. The condition ki = K2 is fulfilled for K2 = VQ = 2.449..., which is indeed where 
Taylor and Herzfeld find the Nb phase. Our choice for K2 implies that our particles have the 
same side ratios close to the Landau point and therefore that the phase-diagram topology 
should be almost identical in that region. This is true, save the discrepancies mentioned 
above. Outside of this region we are exploring particle sizes differently. In fact, in the 
following we go beyond Taylor and Herzfeld's work and consider other values of K2] we will 
see that changes in the phase diagram are very drastic and not just quantitative. 

We have obtained the phase diagram for the case k,2 = 5, which is plotted in Fig. [3]^a). 
New features are apparent. First, the Nb phase again appears for ni ~ K2, but now the 
topology of the phase diagram in this region changes dramatically [a zoom is provided in 
Fig. Et^b)]. Now the Nb becomes disconnected from the I-Nu transition, and there appears a 
direct, first-order transition between the two unaxial nematics, N^^ and N^^ This remark- 
able result, which has never been observed before in the context of molecular theories for 
one-component hard anisotropic particles, resembles the phenomenology obtained by Taylor 
and Herzfeld and discussed above, but with the important difference that their reported 
Nu''-Nu'' transition is continuous (it is likely that, for larger values of K2, their theory would 
have predicted a first-order transition as well). The change in topology of the present case, 
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FIG. 4: (a) Branches of the free energy density per unit thermal energy, pj- /V , as a function 
of packing fraction 1] for the case K2 = 5, ki = 4.76 = 0.21) and for the different phases. 

(p) (r) 

Continuous curve: Nb- Dashed curve: W^'. Dotted curve: NJ)^ Dash-dotted curve: I. The 
symbols indicate bifurcation points, as follows: open square, I-nJp^; rhombus, n{p^-N{j'^ (note that 
the last two transitions are of first order but very weak); open circle: N^^-Nb; filled circle: Nb-S^^; 
filled triangle and square: bifurcation from (unstable) Nu phases to a spatially nonuniform phase, 
(b) Orientational order parameters as a function of packing fraction for the same particle size 
parameters as in (a). The symbols indicate the discontinuity values of Q caused by the first-order 
nature of the Njj^-N^^ transition. 



with respect to the one obtained for the case 1^2 = 2.5, can be explained using a Landau 
expansion for the free energy up to the fourth order in the order parameter tensor Qap- 
Imposing particular conditions on the expansion coefficients, the continuous Nb-NIj''''' tran- 



sitions can be replaced by a first-order N 
on this point). 



u-'-Ny'* transition (see 



for a detailed discussion 



The manner in which the first-order N^'^-Ny'* transition connects to the biaxial nematic 
phase is also peculiar. As the packing fraction is increased following the Ny^'-N^^ binodals, 
first the changes to Nb at a continuous transition, while the transition gap decreases and 
eventually vanishes, see Fig. [^b). Then the transition reappears as a first-order transition 
and again becomes continuous at a tri critical point. These phenomena occur in a very small 
region of the phase diagram. 

The spinodal curves in the case k,2 = 5 present some important changes with respect 
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to the previous case. One is that, due to the higher aspect ratio, the PS phase is very 
unstable. For high ni there appear unstable fluctuations in the N^^ (not the I) phase, and 
these fluctuations are now of columnar symmetry, Cy\ This symmetry involves particles 
arranged with their shortest axes parallel to the columns. As the aspect ratio increases from 
unity, columns become unstable, and the spinodal changes to Sy\ For still larger ki this 
spinodal meets the biaxial nematic region, and reappears in this region as a S^^-' spinodal. 
This phase consists of particles arranged in layers with their shortest axes perpendicular to 
the layers, and with their long axes oriented in the layer planes. Again the spinodal changes 
its nature and joins the spinodal curve of the Sg phase; this is a biaxial smectic phase 
with the long particles axes perpendicular to the layers and with the secondary particle axes 
oriented in the layer planes. Therefore, as the aspect ratio K2 increases, the smectic phase 
tends to develop two types of symmetries, namely rod- and plate-like symmetries. These 



phases were not observed in [3l| due to the low aspect ratio chosen. 
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FIG. 5: (a) Phase diagram in the plane rj — ki for hard biaxial parallelepipeds with K2 = 10. (b) 
A zoom of the region about k^^ ~ 0.1. Labels as in Fig. [2j 

Finally, we note that the region where the Nb phase is stable has considerably increased 
with respect to the K2 = 2.5. Nevertheless, we must bear in mind that the nematic-to-smectic 
transitions could be of flrst order, and therefore that the stability region of the Nb phase 
predicted in Fig. [3t^a) could be much reduced or even vanish. Also, the spinodal lines shown, 
which bound the Nb region from above, have been calculated from the stable phases at the 
particular density chosen. But free-energy differences between unstable uniaxial and stable 
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biaxial nematic phases are very small, as can be seen in Fig. Hl^a). In fact, the unstable N^'^^ 
phases bifurcate [filled triangle and square in Fig. lU^a)] to a non-uniform phase before the 
stable Nb phase does (filled circle), i.e. inside the shaded region labelled Nb in Fig. Hl^a). It 
could happen that the branch to which these unstable Nu'^'' phases bifurcate become more 
stable than the Nb phase before its corresponding bifurcation point. This scenario would 
reduce (but not completely eliminate) the shaded region. For the sake of completeness and 
illustration, the behaviour of the order parameters along the stable free-energy branch in 
plotted in Fig. Ht^b); in particular, we see why A is a more convenient parameter than 6 to 
represent biaxial order in this system. Note the discontinuity in Q indicated by the symbols, 
caused by the first-order nature of the N^^-Ny^ transition. 

Our final result concerns the case K2 = 10, Fig. [5l Here the N^^'-N^^ transition remains 
of first order and again ends in a Landau point as it joins the I-N^'^'' transition. The 
transition gap is considerably decreased, and this is probably related to the discretization 
of particle orientations, which promotes strong orientational ordering of particles with high 
aspect ratios. But an important new feature is that now there appear two biaxial nematic 
phases, Ng^ and Nj^'', coexisting in a narrow density range [see Fig. [S](b)]. The spinodal 
curves appear at more or less similar values of the packing fraction, but the nematic stability 
begins at much lower values, with the consequence that the predicted stability region of the 
biaxial nematic phase is much increased. 



IV. CONCLUSIONS 



Since the motivation of the present work was the observation of biaxial nematic phases 



in goethite 28||, it is natural to ask if our results reproduce the experimental findings and 



if our theory could therefore be used to predict phase behaviour in these systems with any 



confidence. In the experiments reported in 38|, which summarise the work done on goethite 



by the Utrecht group up to now, is in the range 2.2 — 3.0 (probably a bit lower if one takes 
into account the Debye screening lengths of the particles in the solution) and 6.2 ^ ^ 3.1. 
Only in the case ki ~ 3 is a biaxial nematic phase found. This ratio corresponds to 

a phase diagram in between those of Figs. [2] and [3] but, according to our calculations, it 
would not correspond to favourable particle size ratios to easily observe the Nb phase. Our 
prediction is that higher values of both ki and K2 (with ^ ^2) should be explored to find 
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stable biaxial nematic phases in a wider concentration range. However, two aspects require 
further consideration in order to obtain more reliable predictions. One is the calculation 
of absolute free-energy minima of the density functional for the non-uniform phases, which 
could give rise to shifted stability regions in case first-order transitions were found. The other 
is the effect of particle size polydispersity (which in the experiments is in the range 0.25 — 
0.55). This effect may be very important, as polydispersity tends to preclude the occurrence 
of phases with partial positional order (leaving wider Nb stability windows) and induce 
demixing phenomena. These considerations point to the need to incorporate particle size 
polydispersity in the calculations if one is to establish closer contact with the experiments. 
However, from the computational point of view, this task is exceedingly difficult, since there 



exist three, in principle not completely independent, polydispersity parameters [37|. Efforts 
in this direction are now being undertaken in our group. For the moment, we hope that 
our present work represents a step forward in understanding the phase behaviour of biaxial 
particles and can help foster the development of FMT studies in the direction of more 
complex particle shapes jsol, 40 1. 
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